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The story so far... 



In the previous parts of this course, we derived the Boltzmann 
distribution: 

N _ £ i 
~Z 



rij = — e kT. (1) 



We applied this distribution to three example systems: a spin-i 
magnetic solid, a collection of 1-dimensional quantum harmonic 
oscillators, and the Maxwell-Boltzmann gas. We saw how to 
derive a number of properties of these systems, including the 
energy, heat capacity, and entropy as functions of temperature. 
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The story so far... 



We derived the Boltzmann distribution using a number of 
assumptions, including: 

• the states formed by interchanging two particles with 
different energies represent distinct microstates; 

• there is no restriction on the number of particles that can 
occupy a given single-particle state. 

In quantum systems, these assumptions are not necessarily 
satisfied. In fact, we saw already in the case of the "rotational" 
heat capacity of diatomic hydrogen that constraints on the 
wave functions of the atomic nuclei led to observable effects: 
these constraints are not accounted in the above assumptions. 

This suggests that to obtain a more complete and accurate 
description of the thermal properties of materials, we need to 
take a more careful look at the quantum properties of the 
constituent particles. 
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Indistinguishable particles 



We shall begin by re-visiting our assumption that interchanging 
the energies of two particles leads to distinct microstates. In 
other words, instead of assuming that the particles can be 
distinguished from each other, we shall make the assumption 
that the particles are indistinguishable. This is a more realistic 
assumption in many cases, and results in a fundamental change 
in the number of microstates accessible to a system. 

Along the way, we shall find that there are two distinct cases, 
depending on whether single particle states can be occupied by 
at most one, or an unlimited number of particles. 

The final result will be two new distributions for describing 
systems composed of indistinguishable particles. However, we 
will obtain the familiar Boltzmann distribution in the "classical" 
limit. 
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Symmetric and antisymmetric wave functions 



Consider a system with two single-particle states A and B. Let 
Vm(^) be the wave function for a particle in state A, and iPb( x ) 
be the wave function for a particle in state B. Now suppose 
that the system contains two particles. If the interaction 
between the particles is very weak, then we can ignore any 
(small) change in the states A and B. If one particle occupies 
state A, and another particle occupies state B, then an 
approximate solution to the Schrodinger equation for the 
complete system is: 

^(xi, x 2 ) = *Pa( x i)^b( x 2), (2) 

where x\ and x 2 are the coordinates of particle 1 and particle 
2, respectively. 
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Symmetric and antisymmetric wave functions 



If we interchange the two particles, then we form a new wave 
function: 

V^Oi, x 2 ) = ^b( x i)^a( x 2)- (3) 

However, if the particles are indistinguishable, then the 
properties of the system must be unchanged if the particles are 
interchanged. Assuming that the phase of the wave function is 
not observable, then ^ can differ from ^ by (at most) some 
phase factor. 

In other words: 

ip (x 1 ,x 2 ) = e l ip(x 1 ,x 2 ), (4) 

where S is a constant. 
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Symmetric and antisymmetric wave functions 



Thus, the operation of interchanging two particles is 
represented by multiplying the wave function by a phase factor, 
e lS . However, if we interchange the particles twice, we get back 
to the exact situation we started with. This implies that: 

e 216 = 1, (5) 

and hence: 

i& — 



e w = ±1. (6) 



Both cases e lS = +1 and e l6 = — 1 are realised in nature: in the 
first case, we say that the wave function is symmetric with 
respect to the interchange of two particles; in the second case, 
we say that the wave function is antisymmetric. Note that the 
wave function only need be symmetric or antisymmetric if the 
two particles are indistinguishable. 
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Symmetric and antisymmetric wave functions 



To enforce symmetry or antisymmetry, we need to modify the 
wave function (2). A symmetric wave function can be written: 

1 

^(xi, x 2 ) = ^= [^a( x i)^b( x 2) + ^a( x 2)^b( x i)] , ( 7 ) 

and an antisymmetric wave function can be written: 

1 

^Oi, x 2 ) = ^= [iPa( x i)^b( x 2) - ^a( x 2)^b( x i)] • (8) 

The factor l/y/2 is required for normalisation of the wave 
function: 

2 

\^( x li x 2)\ dx\dx 2 = \. (9) 

Note that we assume that the states A and B are orthonormal: 

IVm( x )| 2 dx = 1, / / 0^(x) / 0 j b( x ) dx = 0. (10) 



Statistical Physics 7 Part 4: The Fermi-Dirac Distribution 



Spin and symmetry 



Particles for which the wave function is antisymmetric under 
interchange of any two of the particles are called fermions. 

Particles for which the wave function is symmetric under 
interchange of any two of the particles are called bosons. 



It can be shown, from quantum field theory, that particles with 
half-integer spin are always fermions, whereas particles with 
integer spin are always bosons. 



In this part of the lecture course, we shall consider the 
statistical mechanics of fermions. We shall consider the 
statistical mechanics of bosons in the next part of the lecture 
course. 
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The Pauli exclusion principle 



Regarding fermions, our first observation is that antisymmetry 
of the wave function implies that each single-particle state 
cannot be occupied by the two particles at the same time. If 
we try to put two particles into state A, for example, the wave 
function becomes: 

1 

The wave function vanishes. 



This implies that for fermions, any single-particle state can be 
occupied by at most one particle: this is the Pauli exclusion 
principle. 
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Wave functions for many particles 



With just two particles, it is easy to write down a symmetric or 
antisymmetric wave function. With more particles occupying 
more states, the wave function must be either symmetric or 
antisymmetric with respect to the interchange of any two 
particles. 



The symmetry or antisymmetry can be enforced by writing the 
wave function as the determinant of a matrix, whose 
components are formed from all combinations of particles and 
states. 
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Wave functions for many particles 



For example, with three fermions (with coordinates x\, x 2 and 
£3) occupying three single-particle states (A, B and C), the 
wave function would be: 



^(x 1 ,x 2 ,x 3 ) = 



V3\ 



V>#Oi) 



ip B ( x 2) i>B( x 3) 



(12) 



The corresponding (symmetric) wave function for a system of 
bosons would be constructed in the same way, but each minus 
sign in the evaluation of the determinant must be replaced by a 
plus sign. 
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Distributions for indistinguishable particles 

Recall that we derived the Boltzmann distribution by looking 
for the distribution that contained the largest number of 
microstates. To do this, we needed to write down an expression 
for the number of microstates in any given distribution. 

AH 

t= , — — (13) 

n\\ ri2- ... 



A distribution is specified by giving the number of particles 
occupying each single-particle state. The number of distinct 
microstates corresponding to a given distribution depends on 
whether the particles are distinguishable or indistinguishable; 
and, in the latter case, on whether the particles are fermions or 
bosons. Thus, there will be different distributions in the case of 
indistinguishable particles (Maxwell-Boltzmann distribution); 
fermions (Fermi-Dirac distribution); and bosons (Bose-Einstein 
distribution). 
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Distributions for fermions 



Let us consider states corresponding to energy levels Each 
energy level can be degenerate; i.e. there can be many distinct 
single-particle states with the same energy. Let the number of 
states with energy ^ be and the number of particles with 
that energy be n^. For fermions, it follows from the Pauli 
exclusion principle that the largest number of particles with 
energy si is g it i.e. the maximum value of is g im 
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Distributions for fermions 



The number of distinct ways of filling gi single-particle states 
with rti particles, with at most one particle per single-particle 
state, is: 

U = -yr^ — Vl . (14) 

nil (gi - rii)\ 

(Recall that m objects can be arranged in a sequence in m! 
different ways). 

Each group of g^ states contributes a factor ti to the total 
number of microstates for a given distribution {n^}. Thus, the 
total number of microstates for fermions is: 

«f «"<}) =n „,,,"!. „,), ■ (15) 
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Most probable distribution for fermions 



Using the principle of equal a priori probabilities, the most likely 
distribution in which to find a system of a large number of 
fermions is that with the largest number of microstates. 
Therefore, we need to maximise t F , or equivalently, Intp, with 
respect to small changes in the energy level populations 

From (15), we have: 

Assuming that ^ and are large so that we can use Stirling's 
approximation, we have: 

Intp « Yl [9i ln ^ - n i lnn i - (9i ~ n i) ln (9i ~ n i)]- ( 17 ) 

i 

The change in Intp with respect to changes dni in the 
populations is: 

' 9i ~ n i 



d(ln t F ) « Yl ln " dn *- ( 18 ) 
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Most probable distribution for fermions 



Changes in the populations of the single-particle states are 
subject to the constraints: 



(19) 



and: 

i i 

where TV is the total number of particles in the system, and U 
is the total energy, both of which are constant. 



Combining equation (18) with (19) and (20), with Lagrange 
multipliers a and (3, we have: 



d(ln tp) « Y 



n 



9i ~ n i 



Hi 



+ a + (3e, 



drio = 0 



(21) 
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Most probable distribution for fermions 



Equation (21): 



d(ln t F ) « 



In 



9i ~ n i 



Hi 



+ a + (3e, 



drii = 0, 



must be true for any set of variations dnf, therefore, we can 
write: 



n l - — — | + a + (3si = 0. 

Tin 



(22) 



Solving equation (22), we find that the populations of the 
energy levels in the most probable distribution are: 



rii = 



9-, 



e -a-(3ei _(_ i ' 



(23) 
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The Fermi-Dirac distribution 



Writing: 

e~ a = B, and 0=~, (24) 

where k is Boltzmann's constant, the most probable 
distribution for fermions (23) can be written: 

ni = f ■ (25) 

BekT + 1 



This is the Fermi-Dirac distribution. The constant B can be 
found from the normalisation condition: 

£"i = £— | = N. (26) 

i i BekT + 1 
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The Fermi-Dirac distribution for a dilute system 

At first sight, the Fermi-Dirac distribution looks quite different 
from the Maxwell-Boltzmann distribution, which we derived 
under the assumption that the particles were distinguishable. 
But consider the case that the population of each energy level 
is much less than the number of available states at that energy: 

U% « 1. (27) 



9% 

In this case, we have from (25): 

Be^ + l»l ; (28) 

and therefore: 

BekT + 1 « BekT. (29) 

The Fermi-Dirac distribution can then be approximated by: 

n . ^| e -e\ (30) 

which has the same form as the Maxwell-Boltzmann 
distribution. 
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The chemical potential 



We see that for a "dilute" system, in which the population of 
each energy level is much less than the number of states 
available at that energy level, the Maxwell-Boltzmann 
distribution provides a good approximation to the more 
accurate Fermi-Dirac distribution. 

We expect to see the Fermi-Dirac distribution differ 
significantly from the Maxwell-Boltzmann distribution under 
conditions in which the population of some of the energy levels 
approaches the number of states available at those energy 
levels. This can occur at low temperature, as we now show. 

We begin by defining a parameter n, with units of energy, such 
that B = e ~^l kT . The Fermi-Dirac distribution is then written: 

ni = -^r • ( 31 ) 

e kT + 1 
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The chemical potential 



Now consider the probability that a particular single-particle 
state is occupied by a fermion. Since, at energy e if there are 
particles distributed over ^ single-particle states, the probability 
that a particular state is occupied is: 



Hi 



9% e k T + 1 



(32) 



Now consider the low temperature limit, T — ► 0. We see that in 
this limit, 



Hi 



9i 



= 1 for Si < ji, (33) 



and: 



9i 



= 0 for Si > fi. (34) 



In other words, in the low temperature limit, all states with 
energy below n are occupied, and none of the states with 
energy above \i are occupied. 
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The chemical potential 



The parameter /x is called the chemical potential. In general, 
the chemical potential will be a function of the variables 
specifying the macrostate, including the number of particles, 
and the temperature. 



The chemical potential in the Fermi-Dirac distribution plays a 
similar role to the partition function Z in the Boltzmann 
distribution. It is determined by the normalisation condition; 
depends on the temperature; and parameterises the 
distribution. 



For the Boltzmann distribution, we found that knowledge of 
the partition function was very useful for deriving thermal 
properties of the system. The chemical potential is similarly 
useful in systems of fermions. 
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The chemical potential 



As an example, consider the case that gi is constant. The 
normalisation condition can be written: 

^ r°° ode , 

i 

where there are gde single-particle states in the (small) energy 
range from e to e + de. 



0 e~kT + 1 



Hence, we find that the chemical potential is given by: 

fi(T) = kT\n (eW? _ ij m (3 6 ) 

As we might expect, the low temperature limit is: 

N 

H(T) -> — as T -> 0. (37) 

9 

In the low temperature limit, the N lowest energy states are 
occupied by particles; since the states are uniformly distributed 
in energy, the energy of the highest-energy occupied state is 

N/g. 
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The Fermi temperature 



Now we define the Fermi temperature T F by: 

/i(0) 



T F = 



k 



(38) 



where /i(0) is the value of the chemical potential /i(T) in the 
limit T^O, and k is Boltzmann's constant. /x(0) is sometimes 
called the Fermi energy, and written e F . The Fermi energy is 
related to the Fermi temperature by: 

e F = kT F . (39) 

Let us plot n/g (the fraction of occupied single-particle states) 
as a function of energy for various temperatures; the result is 
shown on the next slide. We see that for: 

T < T F , (40) 

the distribution is significantly different from a 
Maxwell-Boltzmann exponential decay; while for: 

T > Tp, (41) 

the distribution starts to approach a Maxwell-Boltzmann 
distribution. 
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The Fermi temperature 
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The chemical potential 



Note that for e = n, we find: 



1 (42) 



e—a o 

0* e~kT + 1 ^ 



At an energy level equal to the chemical potential, the average 
population of the available states is 1/2. This is a general 
result for the Fermi-Dirac distribution; it is not restricted to the 
case where the density of states is constant. 
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The chemical potential 



It is often useful to know how the chemical potential varies 
with temperature. In the case where the density of states is 
constant, the chemical potential is given by (36): 

N 



H(T) = kT\n (e9 kT - . 



We plot /i(T) as a function of temperature T on the next slide; 
the significant feature is that at low temperatures (T less than 
about Tp/2), the chemical potential is almost constant: it 
varies very little with temperature. 



Note that as the temperature increases, the chemical potential 
Ijl{T) becomes negative. We saw from equation (42) that states 
with energy equal to the chemical potential have an average 
population of 1/2. Thus, when the chemical potential becomes 
negative, all states have an average population less than 1/2. 
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The chemical potential 
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The Fermi gas 



Let us work through an example more thoroughly. Since we are 
going to invest some time and effort, we might as well choose 
an example with some physical significance. An appropriate 
example is provided by the Fermi gas: this is the fermion 
equivalent of the Maxwell-Boltzmann gas. A Fermi gas consists 
of a collection of weakly-interacting fermions that are free to 
move within a bounded volume V, but are unable to move 
beyond the boundary. 



Fermi gases occur in a number of physical situations, including 

• conduction electrons in a metal; 

• liquid 3 He; 

• neutron stars. 
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The Fermi gas: energy distribution 



Let us first find the energy distribution of particles in the Fermi 
gas. At high temperatures, this will approach the energy 
distribution of particles in the Maxwell-Boltzmann gas; but at 
low temperatures, we expect the fermion nature of the particles 
to be important. 



The energy distribution of particles in a Fermi gas is given by: 

n(s) de = g(e) /(e) de, (43) 

where f(e) is the Fermi function'. 

ZOO = ~e=£ ' (44) 

e kT + 1 

and g(e) is the density of states, given in this case by: 

3 

g(e) de = (2s + l)^ (-^) V~e de, (45) 

where s is the spin of the particles in the gas. 
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The Fermi gas: chemical potential 



To find an explicit expression for the energy distribution n(s), 
we need to find an expression for the chemical potential in 
terms of the temperature, T. The chemical potential \i is 
determined by the normalisation: 

roo 

/ n(e) de = N. (46) 

Jo 

Thus, to find an expression for the chemical potential fi, we 
need to solve: 

This is not an easy integral to perform, so we resort to some 
"tricks", as we now proceed to explain. 

To simplify things (slightly) we shall consider the specific case 
of particles with spin s = \. 
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The Fermi gas: chemical potential 



To perform the integral (47), we first consider an integral of 
the form: 

JO de 



Performing this integral by parts gives: 

df(e) 



I = _G(0) - r G(s) de. (49) 

JO de 



Now we use the fact that, for temperatures T <c T Fl the Fermi 
function f(e) is essentially a step function. Thus, the derivative 
vanishes everywhere, except at energies close to the chemical 
potential, where it becomes large: in fact, as the temperature 
tends to zero, the derivative of the Fermi function tends to the 
Dirac delta function. 
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The Fermi gas: chemical potential 



Let us make the assumption that T <cT f . Then, because of 
the properties of df /de, only values of G(e) fove^n contribute 
to the integral /. That means we can expand G(e) in a Taylor 
series around e = n'. 

G(e) = G(/x) + (e - /x)G'(/x) + \{e - fi) 2 G /f ^) + ... (50) 



Substituting this series into equation (49) gives: 

/ = -G(0) + G(/i)/ 0 + kTG f (^)h + ^^G"(/x)J2 + (51) 



where: 
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The Fermi gas: chemical potential 



Now we use: 



e—\i 

df e kT 



^ r eW+l) 2 kT 



(53) 



to give: 



X 



Igtt I x ,,2 (54) 



— oo 



(e x + 1) 



We have again made use of the fact that the integrand is large 
only over a narrow range, this time to extend the lower limit of 
the integration to — oo, and change the variable of integration. 

We can look up the values of Ig for t = 0, 1,2... We find: 

/ = -G(0) + GO) H (fcT) 2 G"(/x) + ... (55) 
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The Fermi gas: chemical potential 



If we choose: 

G(s) = |J, (56) 

so that: 

G'(e) = si, and G"(e) = ^e~i, (57) 
then, applying the result (55) to equation (47), we find: 



3 3 

V /2m\2 r°° xfe V (1m\l (l 3 tt 2 /f ^ 2 -M 



l 



2tt 2 V ft 2 / Jo e ^ + i 2tt 2 V ft 2 J \ 3 12 

(58) 

The approximate solution to equation (58) is: 



1 - 



TT 2 /T^ 2 ' 



12 \T 



(59) 



where kT F = ep, and the Fermi energy e F is given by 



h 2 (3tt 2 nY , an . 
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The Fermi gas: chemical potential 



We can check our result for the Fermi energy (60) by using the 
fact that in the low temperature limit, all single-particle states 
with energy below the Fermi energy e F are filled, and all 
single-particle states with energy above hq are empty. Thus: 



rep 

/ g( £ ) de = N. 
JO 



(61) 



With the density of states g(e) given by (45): 



g(e) de = 



V 



2m\ 2 



yfs de, 



2tt 2 V h 2 

this gives the same result as in (60) above. 



(62) 
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The Fermi gas: chemical potential 



As in the case of a system with constant density of states, we 
find that for T «T f) the chemical potential has a weak 
dependence on temperature (as before, we define kT F = s F ): 
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The Fermi gas: energy distribution 



We have derived explicit expressions (59) and (60) for the 
chemical potential, valid for T <c Tp\ 



/i «3 Ep 



1 - 



7T 



T 



12 \T 



F 



h 2 f3ir 2 N\ 3 



e F = 



2m 



V 



The distribution of particles with energy in the Fermi gas can 
then be written: 



V 



2m\ 2 



2tt 2 \h 2 



£ — fJL 

e kT + 1 



(63) 
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The Fermi gas: energy distribution 
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The Fermi gas: total energy 



The total energy in the Fermi gas is given by: 



roo 

U = en(e) de. 
JO 



Using equation (63) for the particle distribution n(s) 



n(s) de 



V /2m\2 ^Je 



2tt 2 \h 2 



£ — fJL 

e kT + 1 



de, 



we have: 



U = 



V 



2m\ 2 



•oo 



3 
£2 



2tt 2 V h 2 J Jo 



e—\i 

e kT + 1 



de 



(64) 



(65) 



The integral can be performed using the previous results, (48) 
and (55), with: 



G(e) = -si. 
v J 5 



(66) 
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The Fermi gas: total energy 



The result for the total energy of the Fermi gas is 



U = 



V 



2m\ 2 



2tt 2 V h 2 



2 5 7T 2 o 1 

-/I2 + —(kT) 2 fl2 



(67) 



Using equations (59) and (60): 



jji Ki e p 



1 - 



7T 2 T 



12 \ T F 



e F = 



H 2 /3tt 2 AT\3 
2m I V 



we can express the total energy as: 



U = -Ne F 
5 e 



1 + 



5ir 2 f T 



12 XT 



F 



(68) 



where, as usual, kT F = e F , and the above results are valid for 

T <C T F . 
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The Fermi gas: total energy 



We can see whether our result (68) for the total energy in a 
Fermi gas is sensible, by making a rough estimate as follows. 
First, we calculate the total energy in the limit T — > 0, by using 
the fact that in the low temperature limit, all energy levels up 
to the chemical potential e F are occupied, and all energy levels 
above the chemical potential are empty. The total energy in 
the limit T — ► 0 is then: 

V /2m\§ [£f 3 

U ° = ^(^) Jo E2de > (69) 

3 



-el, (70) 



27r 2 \h 2 J 5 
|iV£ F . (71) 



In the last line, we have used, from (60): 



V /2m\§ | , 

^ = ^5 hr) 4- (72) 
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The Fermi gas: total energy 



To estimate the temperature dependence of the total energy of 
the system is a little more difficult. Let us assume that: 

1. at temperature T, particles within energy kT of the Fermi 
energy e F are excited into higher energy levels; 

2. the particles excited into higher energy levels have an 
energy gain |/cT. 



The first assumption is justified by inspection of the plot 
showing the energy distribution of a Fermi gas at different 
temperatures (reproduced on the next slide). The second 
assumption is justified by the principle of equipartition of 
energy, with the further assumption that the particles have 
three degrees of freedom. 
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The Fermi gas: total energy 




r = o.oir F 



s S F 



* ">H„ 





\ 
\ 

* 

* 
♦ 

v 



T = 0.57> 



\ T = 0.1T F 



"•<><„ 



1.2 



1.4 



1.6 



1.8 



Statistical Physics 



44 



Part 4: The Fermi-Dirac Distribution 



The Fermi gas: total energy 



The increase in energy at a temperature T (<c T F ) is therefore 
given approximately by: 



A/7 « g(s F )kT • -kT. 



(73) 



#(£p) is the density of states at the zero-point chemical 
potential e F . For a spin-| Fermi gas, the density of states is 
given by (45): 



Hence, 



which gives: 



, x V /2m\2 I 



3 A/" 




9 / T 
4 iTp 



(74) 



(75) 



(76) 
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The Fermi gas: total energy 



Hence, our estimate for the total energy in the Fermi gas is, for 



3 9 T 

U = U 0 + AU^ -Ns F + -Ne F — 

b 4 \J-F 



(77) 



This agrees quite well with our more careful (but still 
approximate) result (68): 



U = -Ne F 



1 + 



5tH / T 



12 \T 



F 
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The Fermi gas: zero-point energy 



As T — ► 0, the total energy of the Fermi gas becomes: 



3 3 h 2 /3tt 2 

U^ Uo = -Ns F = -— 

5 5 2m V V 



ATI. 



(78) 



This can be compared with the total energy of the Maxwell-Boltzmann gas 
in the same low temperature limit, which we found in a previous part of this 
course: 



C^0,M-B — 3 



h 2 



7T 



2m V V 



N. 



(79) 



The most significant difference between these expressions is that the energy 
of the Fermi gas in the low temperature limit has an additional factor jvi. 
Since N can be large, the zero-temperature energy of the Fermi gas can be 
much larger than the corresponding Maxwell-Boltzmann gas. This is a 
direct consequence of the Pauli exclusion principle, which prevents the 
fermions in the Fermi gas from all falling into the single-particle state with 
the lowest energy. For particles in a Maxwell-Boltzmann gas, there is no 
such restriction. 
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The Fermi gas: heat capacity 



From the temperature dependence of the total energy, we can 
easily find an expression for the heat capacity at constant 
volume: 

Cy = m =^Nk^. (80) 



Again, this result is valid for T <c T F . It shows that at 
temperatures much less than the Fermi temperature, the heat 
capacity is proportional to the temperature, and vanishes in the 
limit T -> 0. 
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The Fermi gas: example of conduction electrons 

We mentioned previously that the conduction electrons in a 
metal may be modelled as a Fermi gas. Let us apply the above 
results to find the Fermi temperature and heat capacity in a 
"typical" case. The Fermi temperature is important, because it 
tells us the range of temperatures for which the expressions we 
have derived are valid. 

The Fermi temperature is given by: 

h 2 /3tt 2 7V\3 

T F = -f-, and e F = ) . (81) 

k 2m \ V 



Note that N/V, the number of electrons per unit volume, is 
given for a monovalent metal by: 

N P N A ^ x 

= ^-A (82) 

where p is the density of the metal, Nj± is Avogadro's constant, 
and A is the relative atomic mass. 
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The Fermi gas: example of conduction electrons 

Let us consider the example of silver, which has one conduction 
electron per atom. Using the following data: 

Density p 10,500 kg m~ 6 

Relative atomic mass A 0.108 kg mol -1 

Reduced Planck's constant h 1.055xl0~ 34 J s 

Electron mass m 9.11xl0 -31 kg 

Boltzmann's constant k 1.38xl0 -23 J K _1 

Avogadro's number 6.023xl0 23 mol -1 



we find: 



e F « 8.82 x 10 -19 J, (83) 



and: 



T F = ^ « 63, 900 K. (84) 

k 

The Fermi temperature is very high by everyday standards; so 
at room temperatures, our expression for the heat capacity 
should work well. 
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The Fermi gas: example of conduction electrons 



At room temperature, T = 293 K <c T F , we can use equation 
(80) for the heat capacity: 

7T 2 T 
Cy = — Nk . 

V 2 T F 



For 1 kg of silver, the number of electrons is 

N = N A /A = 5.58 x 10 24 . Hence, at 293 K, we find that the 

specific heat capacity of the conduction electrons in silver is: 



c v « 1.74 J kg -1 K -1 . (85) 



This compares with the measured specific heat capacity (at 
constant pressure) of silver, 235 J kg -1 K _1 . Clearly, if our 
model is correct, the conduction electrons make only a small 
contribution to the specific heat capacity. 
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The Fermi gas: example of conduction electrons 



Our conclusion is correct: in fact, conduction electrons in 
metals make only a small contribution to the heat capacity of 
the material. We shall consider the dominant contributions, 
from vibrations of atoms in the material, in the next part of 
this course. 

However, as we shall see, the heat capacity associated with 
atomic vibrations varies as the cube of the temperature, 
whereas the heat capacity associated with the conduction 
electrons varies linearly with the temperature: 

C v ,m = 1? + AT 3 , (86) 

where the constant 7 characterises the contribution to the heat 
capacity from the conduction electrons, the constant A (not 
the atomic mass!) characterises the contribution from atomic 
vibrations, and the subscript m denotes the heat capacity per 
mole. 
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The Fermi gas: example of conduction electrons 

This suggests that if we can achieve sufficiently low 
temperatures, the contribution from the conduction electrons 
can start to dominate. 



The heat capacity per mole of the conduction electrons is, 
from equation (80), given by: 

tt 2 ,. . T 

Hence, 



C v ,m = -z-N A k—. (87) 



7 = Mm = l!^. (88) 

T^O dT 2 Tp 

For silver, we find: 

7^ 0.643 mJ mol -1 K" 2 . (89) 

7 is a measurable quantity. We can obtain 7, for example, as 
the intercept of a plot of C v ^ m /T vs T (or T 2 ). The agreement 
between theory and measurement is good, though not exact. 
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The Fermi gas: example of conduction electrons 
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Figure 12.7 — Low temperature molar heat capacity of metals; C Vm /T 
is plotted against T 2 \ The results are consistent with equation (12.29), 
according to which C V/n is equal to yT + /iT 3 . The intercepts give the 
values of y. The values of A can be determined from the slopes. 
(Reproduced with permission from Statistical Mechanics and Properties 
of Matter by E. S. R. GopaL Ellis Horwood Ltd., Chichester.) 
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The Fermi gas: equation of state 



Finally for the Fermi gas, let us find the equation of state, 
which expresses the relationship between pressure, volume and 
temperature. We start from the first law of thermodynamics, in 
the form: 



dU = TdS -pdV. 



(90) 



Thus, the pressure is given by: 



p = 




(91) 
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The Fermi gas: equation of state 



The total energy is given by: 

U = Y j n i e h (92) 

i 

where is the number of particles in a single-particle state 
with energy 

Now, since the number of accessible microstates depends on 
the populations of the states (see equation 15); and the 
entropy (by Boltzmann's equation) depends on the number of 
accessible microstates, the populations are constant at 
constant entropy. Thus: 

^ = Vr3 (93) 

dVJs , l dV' V J 
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The Fermi gas: equation of state 



Recall that, for particles bounded in a box, the single-particle 
energy levels are related to the volume by: 



si oc V 3. (94) 



Hence, 



dsi 2dV 

77 = _ 3^' 



(95) 



and so: 



dV 3V 



(96) 



Thus, summing over energy states, we obtain: 



V = 



9U ) =Y . 2 nA= 2 -- (97) 
Js ,3 V 3 V 



dVJs % 
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The Fermi gas: equation of state 



Using equation (68) for the total energy, we find 



p V = -Ne F 
5 



1 + 



5?H / T 



12 It 



F 



(98) 



This equation of state is valid in the regime T <c T F , and is 
quite different from that for a Maxwell-Boltzmann gas. There 
is a large pressure even in the limit T —> 0; this corresponds to 
the large internal energy that exists in the system, even in the 
limit of low temperature. 



For temperatures much larger than the Fermi temperature, 
T^>T F , the Fermi-Dirac distribution approaches the 
Maxwell-Boltzmann distribution, and the equation of state 
tends towards that of the ideal gas: 

P V = NkT. 



(99) 
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The Fermi gas: equation of state 
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Fermi gas of conduction electrons: magnetic susceptibility 



Finally, we consider the effect of an external magnetic field on 
the conduction electrons in a metal. We again assume that the 
conduction electrons can be modelled as a Fermi gas. Since 
electrons have a magnetic moment: 

Ms = (100) 

2m 

the immediate effect of an external magnetic field is to change 
the energy of the spin states: 

As = ±n B B, (101) 

where the energy of the state increases if the spin associated 
with that state is opposite to the magnetic field, and decreases 
if the spin is parallel to the magnetic field. 
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Fermi gas of conduction electrons: magnetic susceptibility 



For simplicity, let us consider the limit T — > 0. Since, in 
practice, we find that T <c Tp at everyday temperatures, this is 
likely to be a reasonable approximation in commonplace 
situations. 



The reason for taking the limit T — > 0 is that in the low 
temperature limit, we know that energy states up to the Fermi 
energy are occupied, and states above this energy are 
unoccupied. 



The external magnetic field changes increases the energy of 
states with spin opposite to the field, and decreases the energy 
of states with spin parallel to the field. Therefore, to maintain 
full occupancy of states up to the Fermi energy, the number of 
particles with spin opposite to the field must be different from 
the number of particles with spin parallel to the field. 
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Fermi gas of conduction electrons: magnetic susceptibility 



The total magnetic moment per unit volume of the conduction 
electrons in a material is given by: 



(102) 



where n-j- and nj_ are, respectively, the number of electrons per 
unit volume with spin parallel to the external magnetic field, 
and the number of electrons per unit volume with spin opposite 
to the external magnetic field. 



When the external field is applied, a number of electrons with 
spin opposite to the field will "flip" to have spin parallel to the 
field; this maintains a uniform occupation of states up to the 
Fermi energy. 
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Fermi gas of conduction electrons: magnetic susceptibility 

The number of electrons that flip their spins is approximately 
5#(£f) • V-bB, where g(e) de is the number of states in the small 
energy range e to e + de. 




n 




ii 
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Fermi gas of conduction electrons: magnetic susceptibility 



Therefore, in the presence of the magnetic field, there will be a 
difference in the number of electrons with spin parallel to the 
field and spin opposite to the field: 

1 

n T _ n i = 2 x ^T7#( £ f) • VB B - (103) 

Hence, the material acquires a magnetisation (magnetic 
moment per unit volume): 

M = ^g(e F )-fi B B-n B . (104) 

Using the density of states (75): 

3 TV 

g(e F ) = ±— (105) 



the magnetisation is: 



3 TV 9 , 

M = -—-n%B. (106) 
2 Vep 
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Fermi gas of conduction electrons: magnetic susceptibility 



The magnetisation is related to the external magnetic field by 
the magnetic susceptibility % m : 



M = XmH, (107) 



and 



B = nH, (108) 



where /x is the magnetic permeability. Assuming % m <c 1, 

B = M0 (l + Xm)H « n 0 H, (109) 

where no is the permeability of free space. Hence, the 
susceptibility can be written: 

M 3 N 2 

Xm^no— = -^0— V B - (HO) 
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Fermi gas of conduction electrons: magnetic susceptibility 



Let us consider aluminium as an example. Aluminium has three 
electrons in its outer shell; we shall assume that two electrons 
have their spins paired, so only a single electron per atom 
contributes to the magnetic properties of the material. We 
shall further assume that the electrons contributing to the 
magnetic properties can be treated as a Fermi gas. 



We shall use the following properties of aluminium: 



Density p 2,710 kg m - ^ 

Relative atomic mass A 0.0270 kg mol -1 
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Fermi gas of conduction electrons: magnetic susceptibility 

The number of electrons per unit volume is then: 

— = p— = 6.28 x 10 28 rrr 3 , (111) 
and the Fermi energy is: 

t 2m \ V I 



Then, using (110), we find for the susceptibility: 

Xm » 1.1 x 10- 5 . (113) 
This compares with the measured value: 

Xm «2.2 x 10- 5 . (114) 



The agreement is not great - we are out by a factor of two. 
But magnetism is a complicated phenomenon, and at least we 
got the right order of magnitude... 
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Fermi-Dirac distribution: summary 



You should be able to: 

• Explain that the wave function for a system of fermions must be 
antisymmetric under exchange of any two particles; and show that this 
leads to the Pauli exclusion principle. 

• Derive the Fermi-Dirac distribution, as the most probable distribution 
for a system of particles in which a maximum of one particle can 
occupy any single-particle state. 

• Show that, for a dilute system, the Fermi-Dirac distribution can be 
approximated by the Maxwell-Boltzmann distribution. 

• Define, and explain the physical significance of: the chemical potential; 
the Fermi energy; and the Fermi temperature. 

• Derive an approximate expression for the chemical potential of a Fermi 
gas, including the temperature dependence to second order. 

• Sketch the energy distribution for particles in a Fermi gas, showing how 
the distribution changes with temperature (up to the Fermi 
temperature). 

• Derive expressions for the total energy and the heat capacity of a Fermi 
gas. 

• Derive the equation of state for a Fermi gas at temperatures small 
compared to the Fermi temperature. 

• Derive an expression for the magnetic susceptibility of a material whose 
magnetic properties are dominated by a Fermi gas of electrons. 
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